サフ ーー 


22 クト Ao の SAO 


理 に FPGA 内 蔵 の 


本 誌 付属 の FPGA 量 板 の ディ ジタル 信号 処理 へ の 活用 例 を 紹介 ⑯ 二 つ の 動作 モー ド に 対応 
する . ここ で は , 音声 帯域 の スペ クト ラム ・ ア ナラ イザ を 製作 本 機 は , FFT 計算 モー ド と スペ クト ラム ・ ア ナラ イザ ・ 
する . Spartan-3 が 内 蔵 する 乗算 器 を 活用 し て FFT(fast モー ド の 二 つ の 動作 モー ド が あり ます . 
Fourier transform) を 行う . JTAG (Joint Test Action FFT 計算 モー ド は , ホス ト ・ パ ソコ ン か ら 送 っ た デー タ 
Group) イン ター フェ ー ス に よる 通信 機能 を 活用 し , 表示 ・ 制 を FFT 演算 し て ホス ト ・ パ ソコ ン に 送り 返す モー ド で す . 
御 に パソ コン を 利用 し て いる . (編集 部 ) これ は 最低 限 の 構成 の ボー ド で 動作 し ます . 
スペ クト ラム ・ ア ナラ イザ ・ モ ー ド で は , マイ ク で 拾っ 

Spartan-3 に は 高 性 能 な ハー ド ウェ ア 乗 算 器 が 内 蔵 さ れ た 音声 を AD 変換 し , FFT し て ホス ト ・ パ ソコ ン に 送り ま 
て いま す . 付属 FPGA 基板 に 搭載 され て いる XC3S50 に は す . これ は , いわ ゆる 音声 帯域 用 スペ クト ラム ・ ア ナラ イ 
四 つ の 乗算 器 が 内 蔵 さ れ て お り , これ を 使っ て 何 か お も し ザ と し て 動作 し ます . 
ろ い ディ ジタル 信号 処理 が で き な い か な と 考え , 高速 フー 
リエ 変換 FFT : fast Fourier transform) を 用 いた 音声 帯 信 コン デン サ ・ マ イク , ADC, OP ア ンプ な ど で 構 成 
域 ス ペク トラ ム ・ ア ナラ イザ を 作り まし だ 写真 1). 本 機 図 1 に 回 路 図 を 示し ます . コン デン サ ・ マ イク と AD コ 
を 「 フー リエ 次 郎 」 と 名 づけ ます . ン バ ー タ の 間 に OP ア ンズ 米国 National Semiconductor 
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に , 付属 FPGA 基板 を だ れ に で も 気軽 に 楽し ん 
で も ら え る よう に , 外 付 け 回 路 が な く て も 動作 する こと を 
前 提 と し まし た . つま り , 付属 FPGA 基板 に , 電源 と 水晶 
発振 器 と JTAG Joint Test Action Group) ケー ブル を つ 
な い だ だ け の 最低 限 の 構成 で 動作 が 確認 で きる よう に し ま 
す . 操作 は ホス ト ・ パ ソコ ン か ら 行 いま す が , ホス ト ・ パ 
ソコ ン と の 接続 を で きる だ け 簡 便 に する た め , 通信 に は 
FPGA の 書き 込み に も 使わ れる JTAG 信号 を 使い ます . め 
ん どう な 配線 は 必要 あり ませ ん . 
また , 次 の ステ ッ プ と し て , AD コン バー タ と OP ア ンプ , 
そし て 多少 の 抵抗 と コン デン サ を つなぐ こと に より , 音声 帯 
域 の スペ クト ラム ・ ア ナラ イザ を 構成 する こと に し まし た . 写真 1 音声 帯域 スペ クト ラム ・ ア ナラ イザ の 外観 
JTAG イン ター フェ ー ス を 使い , パソ コン で 制御 ・ 表 示す る . 
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TCK 
NJU7223F33 2.5V 


付属 FPGA 基板 較 


水晶 発振 器 隊 


40MHz 較 


表 1 


設計 し た FPGA 機能 の 概要 


内 容 


512 ポ イン ト 


16 ビ ッ ト 固定 
小数 点 


使用 Block RAM 


回 


使用 乗算 器 


固 


FFT 速度 


691.2ws 


ホス ト と の 通信 
方 法 


USER1 命 令 に 
よる JTAG 通信 


スラ イス 使用 率 


44% 
( 768 ス ライ ス 中 


341 ス ライ ス ) 

使用 |/O 数 7 本 

最大 動作 速度 136.129MHz 
任意 形状 の 

その ほか の 

特記 事項 窓 関数 を 処理 


可能 


図 1 音声 帯域 スペ クト ラム ・ ア ナラ イザ の 回 路 図 


社 の LM324) を 入れ , ゲイ ン を 調整 で きる よう に し て いま 
す . A-D コ ン バ ー タ は 米国 Microchip Technology 社 の 
MCP3204 12 ビ ッ ト , 100kHz) を 使い ます . MCP3204 は , 
0V^ 約 4V の 電圧 を , 100kHz の サン プリ ング 周波 数 で 12 
ビッ ト の 2 進 整数 に 変換 し ます . 

写真 1 の よう に , 本 機 は ユニ バー サル 基板 を 使っ て 製作 
し まし た . FPGA ボー ド の ほか に , 3.3V レ ギュ レー タ , OP 
アン プ , A-D コ ン バ ー タ , 電解 コン デン サ な ど を 実装 し て 
いま す . 裏面 に は いく つか の チッ プ 部 品 が 載っ て いま す . 


@ FPGA に 16 ビ ピット 固定 小数 点 FFT を 実装 

本 機 の 中 心 と な る FPGA の 仕様 を 表 1 に 示し ます . 

FPGA では, 512 ポ イン ト の 16 ビ ッ ト 固定 小数 点 の FFT 
を 行い ます . FFT で は バタ フラ イ 演 算 と いう 計算 を 繰り 返 
し 実行 し ます . 512 ポ イン ト の 場合 は , 全部 で 2304 回 の バ 
タフ ライ 演算 を 行う こと に な り ま す . 

この 設計 で は , 1 回 の バタ フラ イ 演算 を 12 ク ロッ ク で 実 
行 す る の で , 1 回 の FFT に は 27648 ク ロッ ク を 要 し ます . ク 
ロッ ク が 40MHz の 場合 は 1 回 の FFT を 約 690 ん s,. すなわち 
FFT だ け な ら は 毎秒 1447 回 計算 で きる こと に な り ま す . 

な お , 詳し い 理由 は 後述 し ます が , デー タ 量 を 512 ポ イ 
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ント と し た の は メモ リ 容量 と の 兼ね 合い に よる も の で す . 
FPGA の アー キテ クチ ャ と デー タ ・ フ ロー を 図 2 に 示し 
ます . バタ フラ イ 演 算 器 と メモ リ ・ ブ ロッ ク を 中 心 と し て , 
全体 を 制御 する ステ ー ト ・ マ シン , FFT 時 に メモ リ な ど を 
管理 する コン トロ ー ラ , ホス ト ・ パ ソコ ン と の 通信 を 行う 


JTAG 通信 回 路 , AD コン バー タ か ら の サン プリ ング ・ デ 
ー タ を 格納 ・ 制 御 す る ユニ ッ ト な どか ら 構 成 さ れ て いま す . 


⑯ パソ コン で 動作 する 制御 ソフ トウ ェ ア を 作成 

パソ コン で 動作 する 制御 ソフ ト ウェア の 画面 を 図 3 に 示 
し ます . Windows 98/2000/XP に 対応 し て いま す . 

FFT 計算 モー ド で は , パソ コン 上 で 作成 し た デー タ を 本 
機 に 送っ て FFT す る こと が で き , サン プル の 関数 を FFT 
し た 結果 を 楽し お こと が で きま す . 

スペ クト ラム ・ アナ ライ ザ ・ モ モー ド で は , FFPT 計算 モー 
ド の 機能 の ほか , マイ ク で 拾っ た 音声 を リア ル タ イ ム に ス 
ペク ト ラム 表示 する こと が で きま す . 


天 FFT 演 算 の く ふ う ) / 


デー タ の 取り 扱い の 利便 性 か ら , 数 値 の 精度 は 16 ビ ッ ト 


付属 FPGA 基板 を 使う た 
回 路 設 計 チ ュ ニ ドリ アル Part2 


16 ビ ッ ト 図 


FFT コ ント ロー ラ 凶 


バタ フラ イ 較 
演算 器 図 


Block RAM 凶 +MX 2 デー タ / 実数 部 図 
未 使用 較 
+/MX 2+1 デー タ / 虚数 部 図 


| 窓 関数 用 図 図 4 複素 数 デー タ の 並び 

窓 関数 デー タ 図 生 算 宮 較 読み 出す メモ リ の アド レス を イン ク 
リ メ ン ト する た びに , 実数 , 虚数 , 
実数 , 虚数 … と デー タ が 出 て くる . 


M 


M 


デー タ 1 虚数 部 
デー タ 2 実数 部 
デー タ 2 虚数 部 


各種 トリ ガ 図 


全体 較 
ステ ー ト ・ マ シン 剛 


図 2 FPGA の アー キテ クチ ャ と デー タ ・ フ ロー 

バタ フラ イ 演算 器 と メモ リ ・ ブ ロッ ク が 機能 の 中 心 で ある . ほか に , 全体 を 
制御 する ステ ー ト ・ マ シン , FFT 時 に メモ リ な ど を 管理 する コン ト ロー ラ , フー 六朗 抽 陣 ブ ログ ラル (CXCopyright 2004 NATIOU RNAi 
ホス ト ・ パ ソコ ン と の 通信 を 行う JTAG 通信 回 路 , Ac-D コ ン バ ー タ か ら の サ ( 

ンプ リン グ ・ デ ー タ を 格納 ・ 制御 する ユニ ッ ト な ど が ある . 


固定 小数 点 と し ます . た だ し , 扱う 数 値 の 範囲 は 絶対 値 が 
1 未満 で ある と し ます . すなわち , Ox0000 を 0 と し , Ox7ff 
を 正 の 数 の 最大 値 09999695), Ox8001 を 負 の 数 の 最小 値 
( - 09999695) と し ます . 土 1 は 使用 し ませ ん . 詳細 は 後述 
し ます が , これ は 計算 の 過程 で 絶対 値 が 1 を 超え な いよ う 


に する た め の 配 慮 で す . 
フー リ エ 変 換 前 の 時 系 列 デー タ も , フー リ 変換 後 の ス 1 _ | II 
ペク ト ル も , 各 点 の 値 は - 1 て 1 の 数 値 に な り ま す . PtarGPXCF が 見 が 9 ま し た 1296 HAD 内 た デー5FFTUC い ます 


図 3 制御 ソフ トウ ェ ア の 画面 
スペ クト ラム ・ ア ナラ イザ ・ モ ー ド の よう す . 


人 @ 複素 数 の 掛け 算 を 行う 
ニニ つ の 複素 数 の 掛け 算 は , 
図 4 の よう な 順序 で メモ リ 上 に 並ん で いま す . 読み 出す メ 
モリ の アド レス を イン クリ メン ト する た びに , 実数 , 虚数 , 
実数 , 虚数 … と デー タ が 出 て きま す . これ は 次 段 の バタ フ 


( gz 太 )X( c 填 の gc- の 十 ( 十 g2)7 
g, ち , c, 9 は 実数 , 7 は 虚数 単位 


と な り ま す . 本 機 で 扱う 数値 の 絶対 値 は 1 未満 な の で , 二 ライ 演算 器 に 入れ る の に つ ご う が 良 く な っ て いま す . 
つの 複素 数 の 掛け 算 を 行っ て も 絶対 値 が 1 を 超え る こと は 
あり ませ ん . gp- pg と pc 十 gg7 に は 加減 算 が 含ま れ て いま @ バタ フラ イ 演 算 と は 
す が , 上 記 の 理由 に よっ て オー バフ ロー する こと は あり ま バタ フラ イ 演 算 は , 現 1) の よう に 三 つ の 複素 数 4, ヵ , 
せん . 間 で 行う 演算 で , 二 つ の 結果 7, が 出力 され ます . 
7 デ 4X gp 


@ 複素 数 デー タ の 並び を 規定 上 


本 回 路 で 使用 する 演算 用 デー タ と 回 転 ひね り 係数 ) は , 


だ 4XWーp ぢ 
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図 5 

バタ フラ イ 演算 器 の 構造 

入力 され る 実数 と 虚数 の デー タ は , 同じ 
バス を 時 分 割 で 共有 し て いる . 


4 メル と ぢ の 絶対 債 は それ ぞ れ 1 未 満 で す が , 加算 と 減算 
の 結果 は 1 を 超え て し まう こと が あり ます . そこ で , 乗算 
と 同時 に 値 を 1/2 に する こと で , バタ フラ イ 演 算 の 結果 の 
絶対 値 が か な ら ず 1 未 満 に 収まる よう に し て いま す . 

これ で は 一 見 , 精度 が 失わ れる か の よう に 思え ます が , 
実は 512 ポ イン ト の FFT で 求め る スペ クト ル は , 最終 的 に 
は 配列 全体 を 1/512 し な けれ ば な り ま せん . FFT 前 に 時 系 
列 に 並ん で いる 各 点 の デー タ は , FFT の 過程 で 計 9 向 の バ 
タフ ライ 演算 を 通る の で , 結果 と し て ちょ うど ぴっ た りな 
の で す 。 


@ 乗算 器 と 加減 算 器 で バタ フラ イ 演 算 器 を 構成 

複素 数 の 乗算 は , 一 般 に は 4 回 の 実数 の 乗算 で 行わ れ ま 
す . 今回 設計 する バタ フラ イ 演算 器 で は , 複素 数 の 実数 部 
分 と 虚数 部 分 を それ ぞ れ 16 ビ ッ ト の レジ スタ に 格納 し て , 
同時 に 計算 で きる よう に し ます . よっ て , 本 回 路 は 2 個 の 
乗算 器 を 使用 し ます . また , 二 つ の 乗算 器 の 後ろ に は 二 つ 
の 加減 算 器 が あり , 乗算 の 結果 を 加減 算 し て 積算 する こと 
が で きま す . 

設計 し た バタ フラ イ 演算 器 の 構造 を 図 5 に 示し ます . デ 
ー タ ・ フ ロー は 図 6 の と お り で す . 入力 され る 実数 と 虚数 
の デー タ は 同じ バス を 時 分 割 で 共有 し て お り , 演算 途中 で 
うま く 処 理 さ れる し くみ に な っ て いま す . 

バタ フラ イ 演算 器 で は 式 1) の 7 と を 求め ます が , 式 
( を 一 見 し た だ け で わか る よう に , 


三 リ アー 2X 


で す . よっ て , 本 回 路 で は いっ た ん J を 求め た 後 , を 2 回 
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MPXIN_WR 
| _ | MPXOUT | 


実数 部 加算 / 減 算 図 一 一 一 
DATAIN 丸め 図 
IPATAIN〉 満 還 MPXIN_AR 、 MR 
の 


ホー ルド 園 LU 虹 


MPXIN_AI 
丸め 凶 


MPXIN_WI 


虚数 部 加算 / 減 算 図 - 


出力 ホー ルド 図 


リセ ッ ト 図 一 ! 


減じ る こと で の 値 を 求め て いま す . 


@⑯ FFT アル ゴリ ズム を 改良 

よく 目 に する 了 W. Cooley と 了 . W.Tukey の FFT アルゴ 
リズ ム は , 図 7 の よう な デー タ の 流れ に な り ま す . これ を 
ソフ トウ ェ ア で 記述 すれ ば すぐ に わか り ま す が , この 演算 

普通 に 行う と 3 重 の for 文 の ルー プ に な っ て し まい ます . 
し か も , 内 側 の 二 つ の ルー プ は 単純 な for 文 ルー プ で は な 
く , ルー プ 回 数 が 毎回 変わ る めん どう な も の に な り ま す . 

この 3 重 の ルー プ を その まま 回 路 で 実現 し よう と する と 
非常 に や っ か いで す . そこ で , バタ フラ イ 演算 の 計算 順序 
と メモ リ へ の 格納 方 法 を 入れ 替え る こと で , 3 重 の ルー プ 
を 1 重 の ルー プ に ほぐす こと を 考え まし た . 改良 し た FEFT 
の アル ゴリ ズム は , 図 8 の よう な デー タ の 流れ に な り ま す . 

この アル ゴリ ズム で は , バタ フラ イ 演 算 の 順序 を 変更 し 
て いる だ け な の で , 演算 回 数 は 通常 の FFT と 変わ ら ず , 
が X log WV/2 の まま で す . し か し , ルー プ が 1 重 に な っ た こ 
と と , 次 に 計算 する べき アド レス や 格納 すべ き ア ドレ ス が 
科 単 に 求まる こと に より , ハー ド ウェア と し て 実装 し や 
まがり ザ 。 

この アル ゴリ ズム の 動作 を 図 8 の よう な 8 ポイ ント FFT 
の 場合 を 例 に 考え て みる と , アド レス 0 番 と 1 番 の デー タ を 
演算 し て その 結果 を 0 番 と 4 番 に 格納 し ます . 次 に 2 番 と 3 
番 の デー タ を 演算 し て その 結果 を 1 番 と 5 番 に 格納 し ます . 
この よう に , アド レス ・ カ ウン タ の ビッ ト ・ シ フト と 最上 
位 ビ ッ ト を ′ 1 に する が 0O' に する か と いう 非常 に 単純 な 
規則 で , 取り 出す メモ リ の アド レス と 書き 込む な メモ リ の ア 
ドレ ス が 求め られ ます. 


% 
DATAD 較 
WBUS 図 
MPXIN AR 
MPXIN_A 較 
MPXIN_WR 
MPXIN_W 
MPXOUT_R 
MPXOUT_ 


ADDOUT_R 


ADDOUT_ 


回 に 図 に 図 に 図 図 に 図 ^ 


た 実数 虚数 図 | 
,) | A: 入力 A : 入力 B 図 
DOUT | W: 入力 W Ki 結果 K 図 


ADDOUTD 


ナ : 結果 」 
本 文 式 1) 参照 較 


図 6 バタ フラ イ 演 算 器 の デー タ ・ フ ロー 
入力 され る 実数 と 虚数 の デー タ が 演算 途中 で うま く 処理 され る . 


また , 回 転 子 の 並ぶ 順 痛 も すっ きり と し て いま す . 例え 
ば , 通常 の 16 ポ イン ト の FEFT の 場合 を 考え て みる と , 
転 子 の 添え 字 ば 00.000000【040 
40 40 4 れれ 02460246 ル 012 
3 4 5, 6, 7}) と な り , 規則 性 を 求め る の は ちょ っ と め 
ん どう で す . 

一 方 , 改良 され た アル ゴリ ズム で 0 0 0 0 0 
00 0 ルル 00004444 4 人 れ 00224 
4 6 6【 0 1. 2 3 4 5 6 7) と いう 順序 に な り 
ます . つま り , アド レス ・ カウンタ の 値 を ビッ ト ・ マ スク 
する だ け で 求め る こと が で き , ハー ド ウェ ア に よる 実装 が 
と て も 簡単 に な り ま す . 

た だ し , この アル ゴリ ズム に は 一 つ だ け 欠 点 が あり ます . 
通常 の FFT で は , 0 番地 の デー タ と 1 番地 の デー タ を 演算 
し た 結果 は , 0 番地 と 1 番地 に 書き 戻し ます . と ころ が , 本 
アル ゴリ ズム で は 読み 出し た アド レス と 書き 込む アド レス 
は 一 致し な い の で , 2 バン ク 以上 の メモ リ を イン タリ ー ブ 
し て 使う こと に な り ま す . また , 2 僅 の メモ リ 容量 が 必要 
と な り ま す . 


回 


111 1 


ステ ー ジ 1 ステ ー ジ 2 ステ ー ジ 3 
図 7 通常 の FFT の デー タ ・ フ ロー 


8 ポイ ント FFT の 場合 . 


ステ ー ジ 1 ステ ー ジ 2 ステ ー ジ 3 
図 8 設計 し た FFT の デー タ ・ フ ロー 
8 ポイ ント FFT の 場合 . 


3. FPGA 内 疫 機 能 の 用 移 


@⑯ BSCAN を 使っ て JTAG 通信 

Spartan-3 に は , JTAG の 信号 を 使っ て FPGA 内 部 の ユ 
ー ザ 回 路 と 通信 で きる よう に する た め , BSCAN コ ン ポ ー 
ネン ト が 用 意 さ れ て いま す . USER1 と USER2 と いう JTAG 
命令 を 使い , TDI, TDO 信 号 を ユー ザ の 用 意 し た デー タ ・ 
レジ スタ に 接続 で きる よう に な り ま す . 

BSCAN の 使用 例 を リス ト 1 に 示し ます . な お , この 例 
で は 多少 簡略 化し て お り , 実際 の 本 機 の 構成 は も う 少 し 複 
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雑 で す . 

BSCAN に は , USER1 と USER2 の 命令 に 対応 し て 二 つ 
の 系 統 が あり ます . BSCAN を 使う た め に は , 任意 の 長 さ 
の シフ ト ・ レ ジス タ を 用 意 し , シフ ト ・ レ ジス タ の 入力 に 
TDI を , シフ ト ・ レ ジス タ の 最 下 位 ビ ピッ ト に TDO を つ な 
ぎ ま す . 

USER1 を 使う 場合 ,。 シ フト ・ レ ジス タ は DRCK1 で シフ 
ト さ れる よう に し ます . FPGA か ら パ ソコ ン に 送信 する デ 
ー タ は , JTAG_ CAPTURE が アサ ー ト され て いる と き の 
DRCK の 立ち 上 が り で シフ ト ・ レ ジス タ に 代入 し ます . ま 


リス ト 1 


BSCAN の 使用 例 


た , UPDATE 信 号 は シフ ト 動作 が 完了 し て デー タ が 確定 
し た こと を 示す 信号 で す . UPDATE の と き に は DRCK は 
出 な い の で ご 注意 くだ さい . 

また , JTAG の クロ ッ ク と シス テム ・ ク ロッ ク は 同期 し 
て いな いこ と に 注意 し て くだ さい . リス ト 1 の 例 で は 
UPDATE 信 号 を シス テム ・ ク ロッ ク で 同期 化し て いま す 
が , 本 来 は UPDATE 信 号 を クロ ッ ク 入 力 と し て 使っ て シ 
フト ・ レ ジス タ か ら 内 部 の デー タ ・ レ ジス タ ヘ デー タ を 移 
す の が 正 し い の で し ょ う . さら に , CAPTURE 信 号 も シフ 
ト ・ レ ジス タ の 非同期 セッ ト に 用 いる の が 正しい の か も し 


ggna1 CLK : gtd 1ogic:-- シ ステ ム ・ ク ロッ ク 

81gna1 SHTFTREG : gtd 1ogio_vector(23 downto 0) : -- シ フト ・ レ ジス タ 
g1gna1 COMMAND TO FPGA : std 1ogioc vector(23 downto 0) : --FPGA に 与え る コマ ンド 
gigna1 STATUS OF FPGA  : std 1ogic vector(23 downto 0) : -- ホ スト へ 送る ステ ー タ ス 
gtgna]1 UPDATE : std 1ogio: -- BSCAN コ ン ポ ー ネ ン ト の uppArg 信 号 

gtqna1 UPDATED : gtd_1ogio,: -- UPDATEg 信 号 を CrK で サン プリ ング し た も の 
stgna1 UPDATED2 : gtd 1ogicz: -- UPDATED 信 号 を CrK で 遅延 させ た も の 


componen モ t BSCAN SPARTAN3 
port ( 


さ 


procesg(CLK) begin 


E (DRCK1'′ evenE and (DRCK1= 1 )) then 


CAPTURE : out STD ULOGTC, ュ 
DRCK1  : out 8TD ULOGTC: ーー シフ ト 動作 は DRCK で 行う 

DRCK2 : out STD ULOGTC: 1E((SEL1 = 1 ) and (CAPTURE = 1 )) then 
RESET 。 : out STD ULOG1C --FPGA 内 の ステ ー タ ス 情 報 を シフ ト ・ レ ジス タ に 代入 

0 out STD OLOGTC。 SHTETREG <= STATUS OF FPGA: 

SEr2 : out 8TD ULOGTC, e1g1f(SEh1 = "1 ) then 

SHTgmT 。 out STD ULOGTC: --SEr1 が アサ ー ト され て いれ ば シフ ト 

敵 。 out STD OLOGTC。 BHTFTREG <= TDT & SHTFTREG(23 downto 1) 
UPDATE  : out STD ULOGTC: end if: 

TDO1 : 3n STD ULOGTC end 1f: 

TDO2 : tn STD ULOGTC 


E(CLK event and (CLK= 1 )) then 
UPDATED <= UPBDATE and SETL1 : 
UPDATED2 <= UPDATED : 
-- シ ステ ム ・ ク ロッ ク で 同期 化す る 


) : 


end componen : 


E ( (UPDATED2 ^ 0′) and (UPDATED = "1 


COMMAND TO_FPGA <= SHTFTREG: 


) ) then 
end FE: 
end 1F: 
end proOoeS8: 
TDO1 <= SHTFTREG(0): 


Block RAM Block RAM 者 Block RAM #2 Block RAM お 
0x0000 - 1 0x0000 
演算 用 デー タ 凶 sin 
に お ける デー タ の 移動 了 ー ュ 
0 _。 cos 
回 転 子 較 
計算 し た い 鐘 | | 計算 結果 の 格納 図 未 使用 図 (sm eos) 図 ooooo 
デー タ の 格納 図 1 
図 9 
Block RAM の 割り 当て 
ー つ の Block RAM に は 512 ポ イン 
が / 


ニニ つの Block RAM を 演算 用 に 使 
う . 窓 関 数 と 回 転 子 を 合わ せ て , 


ー つ の Block RAM を 使用 する . 偶数 回 目 の ス テー ジ に お ける デー タ の 移動 図 
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れ ま せん . 

JTAG の し くみ や , JTAG の デー タ ・ レ ジス タ と は 何 か 
と いう こと に つい て は , 参考 文献 5) を 参照 し て くだ さい . 
な お , BSCAN の 概要 に つい て は , ISE WebPACK が イン 
スト ー ル され て いる ディ レク トリ の 中 の ぎ doc\japanese\ 
de\libs に bscan_spartan3pdf と いう ファ イル が あり ます . 


人 @ 三 つ の メモ リ ・ ブ ロッ ク (Block RAM) を 利用 

XC3S50 に は 四 つ の メモ リ ・ ブ ロッ パ Block RAM) が 内 
蔵 さ れ て いま す . 一 つの サイ ズ は 18K ビ ッ ト な の で , 18 ビ 
ッ ト 幅 で 1K ワ ー ド の メモ リ と し て 使用 する こと が で きま す . 

本 設計 で 扱う デー タ は 16 ビ ッ ト 幅 で す . 一 つの Block 

に は 1.024 ワ ー ド の デー タ を 入れ る こと が で きま す が , 
複素 数 で 計算 する の で 一 つの デー タ に つき 2 ワー ド 使い ま 
す . よっ て , 一 つの Block RAM に は 512 ポ イン ト の 複素 数 
デー タ が 格納 で きる こと に な り ま す . 本 機 で は , 改良 され 
た FFT ア ル ゴ リ ズム を 用 いる た め , Block RAM は 2 バン 
ク 必 要 に な り ま す . し た が っ て , 二 つ の Block RAM を 演 
算 用 に 使い ます . 

回 転 子 は Wo …, Wsss の 256 個 の 複素 数 な の で , 全部 
で 512 ワ ー ド を 使い ます . 窓 関数 は 512 ポ イン ト 分 あり ます 
が , こち ら は 実数 値 な の で 512 ワ ソー ド で すみ ます . 窓 関 数 
と 回 転 子 を 合わ せ て 16 ビ ッ ト 1024 ワ ー ド に な る の で , バ 
ンク 3 の Block RAM に ぴっ た り 納まり ます . 

以上 を まとめ る と , Block RAM の 割り 当て は 図 9 の よう 


0) FPGA の デバ ッ ク グ 法 


オシ ロス コー プ や ロジ ッ ク ・ ア ナラ イザ が な いと FPGA 回 路 
の デバ ッ グ は と て も た い へ ん で す . JTAG ケー ブル で FPGA の 端 
子 の 状態 を 調べ た り , ロジ ッ ク ・ ア ナラ イザ の よう に 波形 を 観 
測 で きる と 便利 で す . ケー ブル を つなぎ か える こと な く , FPGA 
や コン フィ グレ ーション ROM へ の 書き 込み も 行い た いも の で す . 
筆者 は , これ ら の 機能 を 持つ FPGA 開発 ツー ル を 開発 し て い 
す . 評価 版 は 下記 の アド レス か ら ダ ウン ロー ド で きま す . 


http://Www.nahitech.com/jtag/ 


図 A 

MITOUJ TAG の 画面 

FPGA の 余っ た ピン に 内 部 信号 を 出し て デパ ッ グ 
し て いる よう す . 


に な り ま す . 


また , Block RAM は 真 の デュ アル ポー ト ・ メ モリ です . 
本 設計 で は ポー ト A を BSCAN コ ン ポ ー ネ ン ト か ら 操 作 さ 
れる ポート と し , ポート B を FFT の 演算 で 用 いま す . この 


乗数 18 ビ ッ ト 較 被 乗 数 18 ビ ッ ト 較 


A 


符号 較 ]。。 ト 図 


回 四 因 田 回 国 回 回 ME 回 回 回 回 回 
ペー - トー ンーv ニャ ルー デジ 
AXE 16 ビ ッ ト ) 図 
( a) 乗算 器 を 用 いた 固定 小数 点数 AXB の 計算 方 法 較 


乗数 18 ビ ッ ト 図 被 乗 数 18 ビ ッ ト 図 
A 


符号 図 本 ト 凶 


回 因 田 四 思 四 回 EEME 加 回 回 回 回 

ペー ニーーvv ココ ヤヤ ゼロ ャ / 
AXB+2 16 ビ ッ ト ) 図 

( b) 乗算 器 を 用 いた 固定 小数 点数 AX B+ 2 の 計算 方 法 凶 

図 10 固定 小数 点 の 乗算 

18 ビ ッ ト どう し の 乗算 結果 か ら 16 ビ ッ ト の 演算 結果 を 得る . 


TATE OP<1 
G5TATE OF の 
BSTATE of や 
BSTATE OPG> 
BSTATE of< み の 
BSTATE_OP<3> 
BFUY_MEMOWROP 
BFLY.MEMIWROP 


BFLYJMEMOOOUT.OR 
LEDGREENOP 
LEORED.OP 


OO MMMMMIMMMMM MMMMIU 昌 MIM 咽 MI 卓 章介 


ADCCtk 
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よう に 分 ける こと で , デー タ ・ バ ス や アド レス ・ バ ス に マ 果 を 2 で 割っ た 数 が ほし いと き に は 図 10 b) の よう に する 

ル チ プ レク サ や アー ビタ を 用 いる こと な く , 両方 か ら 独立 どう まく の く よ ラ で す 。 

に アク セス する こと が で きま す . これ を VHDL で 表す と リス ト 2 の よう に な り ま す . 
Block RAM2 は 本 機 で は 使用 し ませ ん . 機能 の 拡張 に ご な お , 実際 の 設計 で は , バタ フラ イ 演 算 器 内 で リス ト 3 

利用 く ださい. の よう に , 切り 捨て られ て る 部 分 の 最上 位 ビッ ト を 見 て 科 


単 な 丸 め 処 理 を 行っ て いま す . 


@ 18 ビ ピット の 乗算 器 で 16 ビッ ト 乗 算 を 実行 


Spartan-3 に 内 蔵 さ れ て いる 乗算 器 は 18 ビ ッ ト 符号 付き 2. ホス ト ・ パ ソコ ン と の 通信 移 
の 乗算 を 行う こと が で き , 入力 され た デー タ は つね に 符号 


付き 整数 と し て 解釈 され ます . ホス ト ・ パ ソコ ン と 本 機 の 通信 を 行う た め に , 24 ビ ッ ト 
16 ビ ッ ト の 符号 付き 整数 の 乗算 を 行う た め に は , 下位 2 の コマ ンド ・ レジ スタ と ステ ー タ ス ・ レジスタ を 用 意 し て 
ビッ ト に 何ら か の デー タ を 詰め なけ れ ば な り ま せん が , 本 いま す . コマ ンド ・ レ ジス タ に 値 を 読み 書き する こと で , 
回 路 で は 最 下 位 ビ ッ ト を コピ ー し て 詰め て いま す . メモ リ へ の アク セス や FFT 開始 な どの コマ ンド を 与え る こ 
18 ビ ッ ト どう し の 乗算 結果 は 36 ビ ッ ト に な り ま す が , 本 と が で きま す . また , ステ ー タ ス ・ レ ジス タ の 値 を 読み 込 
回 路 で 欲し い の は 16 ビ ッ ト の 演算 結果 な の で , 下位 の 何 ビ むこ と で , メモ リ の 内 容 や FEFT の 実行 状況 を 知る こと が で 
ッ ト か は 切り 捨て る こと に な り ま す . きま す . な お , コマ ンド ・ レ ジス タ と ステ ー タ ス ・ レ ジス 
本 回 路 で 使う 固定 小数 点 の デー タ 形式 で は , 図 16 a) の タ は JTAG 経由 で シリ アル に アク セス され る の で , コマ ン 
よう に する こと で 乗算 が で きる よう で す . また , 乗算 の 結 ド 書き 込み と ステ ー タ ス 読 み 出し は 同時 に 行わ れ ま す . 


@ コマ ンド ・ レ ジス タ と ステ ー タ ス ・ レ ジス タ 


リス ト 2 18 ビ ッ ト ど うし の 乗算 結果 か ら 16 ビ ッ ト の 演算 結果 を 得る コマ ンド の 形式 を 図 11 に , ステ ー タ ス の 形式 を 図 12 に 
VHDL コー ド 示し ます . 

図 13 は , JTAG 経由 の メモ リ ・ ア クセ ス 回 路 の 原理 的 な 
gae BN gate seaclhs SRCaRS 3 構造 で す . FPGA 内 に アド レス ・ カ ウン タ が あり , メモ リ 
3532R64 80208C8CH5IGKS8RBgEIHS382wmegs014 3 の 読み 書き は この カウ ンタ の 示す アド レス に 対し て 行わ れ 
port ( ah dm 5HD 2ooFC Wacro (0820 ます . カウ ンタ は メモ リ の 読み 書き に よっ て 自動 的 に イン 
RNN78R2R188:929RKak クリ メン ト され る ほか , コマ ンド の ビッ ト 20 を 1T に する 

end component: と 値 を セッ ト で きま す . 
読み 出さ れ た デー タ の ビッ ト 15~ て ビッ ト 0 は , メモ リ の 
GAAS 値 か , も し く は AD コン バー タ か ら 得 らち れ た 現在 の アナ ロ 
RM グ 値 0 て 4095) で す . A-D コ ン バ ー タ の 値 を 読み 込む に は , 
LBPRED ま HPRT3I な あら か じ め ビ ッ ト 22 を "1 に し た デー タ を コマ ンド ・ レ ジ 
SPP ご に 3 スタ に セッ ト し て お きま す . ビッ ト 22 が O の 場合 は メモ 


DATA P2 <= P(35 downto 20): -- AxB/2 の 値 


リ の 内 容 が 読み 出さ れ ま す . 


リス ト 3 -- AxB/2 を 求め る 
丸め 処理 1E (MPXOUT R P(19) = ′1' ) then -- 乗 算 結果 の 19 ビ ッ ト 目 を 見 る 
-- 繰 り 上 が り 
MPXOUT _R <= MPXOUT R P(35 downto 20) + "0000000000000001”: 


e1ge 
-- 繰 り 上 が りな し 
MPXOUT _R <= MPXOUT R P(35 downto 20) 
end 1F: 
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参考 まで に , Block RAM の 内 容 を 読み 出す コー ド を 
リス ト 4 に 示し ます . 
FFT を 実行 する に は , コマ ンド ・ レ ジス タ の ビッ ト 23 を 


に し ます . この 場合 , FFT の 実行 に は 二 つ の モー ド が 
あり ます . 一 つ は ホス ト ・ パ ソコ ン か ら 時 系 列 デ ー タ を 
Block RAM に あら か じ め 書 き 込ん で お く 方 法 , も う 一 つ 
は A-D コ ン バ ー タ で サン プリ ング され た デー タ を FFT する 


MSB 


回 回 回 回 回 四 四 回 回 回 回 四 還 EEEEEEIEEEHE 


方 法 で す . 後者 の サン プリ ング ・ デ ー タ の 場合 は , ビッ ト 
21 と ビッ ト 23 を "1 に し て コマ ンド を 与え ます . ビッ ト 23 
を 1 に し た と き の ビ ッ ト 6 て ビット 0 は サン プリ ング 周期 
と な り ま す . この 値 を 了 了 と する と , サン プリ ング 周波 数 は 
100kHz( 7 十 1) と な り ま す . 

な お , 読み 出し た ステ ー タ ス の ビッ ト 23 が 1 の 場合 は , 
現在 , FFT の 処理 を 行っ て いる 最 中 で ある こと を 示し ます . 


LSB 


L 


書き 込み デー タ ま た は メモ リ ・ ア ドレ ス ・ サ ンプ リン グ 周 期 較 
メモ リ ・ バ ンク 選択 ' 00" 衣 バン ク 0,“" 01" 製 バ ンク 1, "11" バ ンク 3 
デー タ を 読み 出し て , アド レス ・ カ ウン タ を +1 す る 図 
デー タ を 書き 込ん で , アド レス ・ カ ウン タ を +1 す る 図 
ビッ ト 9 て 0 の 値 を アド レス ・ カ ウン タ に セッ ト する 図 


テ A- じ 変換 の 値 を サン プリ ング し て , FFT を 行う 較 


図 11 


テ 1 の 場合 , メ モリ の 内 容 で は な く , 現在 の A-D 変 換 の 値 を 取得 する 凶 


コマ ンド ・ レ ジス タ 


MSB 


= FFT 開 始 、 リ が サン プリ ング 周期 7 を ビット 6 ビッ ト 0 に 設定 する ) 凶 
サン プリ ング 周期 数 =100kHz/(7+1) と な る 較 


LSB 


回 回 回 回 回 四 田 四 加 回 回 四 還 EEEEEEIEEEHE 


( 固 


図 12 


01” 図 読み 出さ れ た メモ リ の 値 , ある い は A- じ 変換 の 値 凶 
定 ) 較 


ステ ー タ ス ・ レ ジス タ 


MEMWR 
ェ ー テ MEMRD 


ー FFT 状 態 Q 待機 中 , + 処理 中 ) 図 


バン ク 選 択 凶 


ADCGET 


書き 込み デー タ 攻 


FFT 回 路 へ 


コン ポー ネン ト 隊 
TDI 


CAPTURE 


MSB 
図 13 0 
JTAG に よる メモ リ ・ ア 

クセ ス 回 路 の 構造 

FPGA 内 に アド レス ・ カ ウ 

ンタ が あり , メモ リ の 読み 


ステ ー タ ス 較 


ADCGET 
Zm_ 9 


書き は この カウ ンタ の 示す 


アド レス に 対し て 行わ れる . 「「「 処 理 中 較 


A-D コ ン バ ー タ 図 
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リス ト 4 
ミ 1d read mem(1inE bank , 1n addgr,1n 1]en, FE 上 1n ヒ 16 *daa) 
Block RAM の 内 容 を 読み Au ー 
内 容 を 読み 出す { 


unsigned char di [25] : 
unsigned char doocap [25] : 
unsigned 1ong ndata: 


//U8ER1 イ ンス トラ クシ ョ ン の 送信 

memset (di ,1,1r 1en) , // すべ て 1 に セッ ト 

tag_setbsreg_From st て (tdi,” 000010”,6) ,// USER1 (000010) を Spartan3 に セッ ト 
tag_sh1ft 1 て (di,tdoocaDp,1 エ 1en) : 


// デ ー タ の 送信 
ndata = Ox00100000 | ((bank & 3) << 16) | (addr & Ox3fFF) : // set memory addressg 


tag_setbsreg_fFrom u1ong (tdi , indata, bsreg_1en, 0): 
tag_sh1ft dr (di ,tdooap,bsreg_1en) : 


// デ ー タ の 受信 

ndata = Ox00040000: // memory read command 
tag_setbsreg_From u1ong (tdi , indata, bsreg _1en, 0): 
For (in 1=0:1<1en : ユ ++) 


{ 


tag_sh1ft dr (di ,tdooap,bsreg_1en) : 
data [1] = ] て ag _getu1ong_ From bsreg (tdoocap, 16, 0) : 


) 


人 @ JTAG で 通信 する ソフ トウ ェ ア 

本 機 は , 米国 Xilimx 社 の Parallel Cable IV を 介し て 
FPGA と 通信 し ます . これ ら の JTAG ケー ブル は パラ レ 
ル ・ ポ ー ト と JTAG 信号 線 が 1 対 1 で つなが っ て いる 構造 な 
の で , 単純 に パラ レル ・ ポ ー ト に 入出 力 す る だ け で アク セ 
AG きま す : 

Windows 2000/XP で 動作 する ソフ ト ウェ ア か ら 1/O ア 
クセ ス を 行う た め に は , 何ら か の デバ イス ・ ド ライ バ が 必 
要 に な り ま す . 今回 は , 参考 文献 6) の PCI デ バッ グラ イ 
ブラ リ for Win32」 を 使用 し に ます. この ライ ブラ リ は , シ 
ステ ム ・ デ ィ レ クト リ ま た は カレ ント ・ デ ィ レ クト リ に 
PCIDEBUGDLL, PCIDBG95.VXD, PCIDBGNT.SYS 
の 三 つ の ファ イル が あれ ば , _IoReadChar ま た は _ 
IoWriteChar を 呼ぶ だ け で I/O ア クセ ス が 行え る も の で す . 
この ライ ブラ リ は Web か ら ダ ウン ロー ド する こと も で き 
まず http://www.ottoto/ kasiwano/newpage16htm). 

この ライ ブラ リ の イン スト 一 ル は 非常 に 簡単 で す が , I/O 
アク セス の た びに カー ネル ・ モ ー ド へ の 切り 替え が 発生 す 
る の で 速度 は 遅く な り ま す . そこ で GIVEIO と いう モジ ュ 
ー ル を 使っ て アク セス する 方 法 に も 対応 し まし た . GIVEIO 
が イン スト 一 ル さ れ て いる パソ コン で は , GIVEIO が 優先 
的 に 使わ れ ま す . GIVEIO は イン スト ーー ル は めん どう で す 
が , カー ネル ・ モ ー ド へ の 切り 替え が な い 分 だ け 高 速 で す . 

本 ソフ ト ウェ ア を 起動 する と , JTAG ケー ブル の 先 に 接 
続 さ れ て いる デバ イス の 構成 を 確認 し ます . デバ イス が 確 
認 さ れる と Block RAM に , 回 転 子 と 窓 関 数 の デー タ を 


間 
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き 込み , ユー ザ か ら の 操作 を 待ち ます . 

スペ アナ ・ モード では, FFT の 後 , 512 ワ ー ド の 演算 結 
果 複素 数 で 256 ポ イン ト ) の 内 容 を 読み 出し ます . これ は 
演算 結果 の 半分 の 量 で す . 一 般 に 実数 デー タ を FFT す る と 
半分 の 周波 数 の と ころ で 折り 返し た 形 の スペ クト ル に な る 
た め , 半分 だ け 読め ば 十分 で す . 

FFT 計算 モー ド で は , ホス ト ・ パ ソコ ン で デー タ を 計算 
し て それ を FPGA 内 の メモ リ に 転送 し , FFT を 行っ て 表示 
まま 

普通 に サン プリ ング し た デー タ を その まま FEFT する と い 
ろ い ろ な 問題 が 生じ る の で , 窓 関 数 と いう 処理 を 施す の が 
一 般 的 で す . 本 機 は 窓 関数 用 の 係数 メモ リ を 持っ て いて , 
この メモ リ に 任意 の 形 の 窓 を 書き 込む こと が で きま す . 本 
ソフ トウ ェ ア は , ハニ ング 窓 , ハミ ング 窓 , 方 形 窓 な ど に 
対応 し て いま す . 


5. 製作 し た 回 路 の 動作 と 性 能 評価 移 


FFT し た デー タ を パソ コン に 取り 込む と , 画面 の 更新 は 
毎秒 8 回 程度 で す . これ は , JTAG 通信 の 速度 が ボ ト ルネ 
ッ ク と な っ て いま す . 

パラ レル ・ ポ ー ト か ら TCK を 操作 する の に 上げ下げ で 2 
ステ ッ プ , TDI の 操作 と TDO の 入力 も 含め る と 4 ステ ッ プ 
を 要 し ます . パラ レル ・ ポ ー ト の アク セス 速度 は , 通常 の 
PC/AT 互換 機 な ら 約 1MHz な の で , デー タ の 転送 速度 は 
250kbps 程 度 た と なり ます . 本 機 で は 24 ビ ッ ト の レジ スタ に 


デー タ を 入出 力 し ます が , TAP の 制御 な ども 含め る と 7 
回 /s5 の デー タ の 入出 力 が 限度 と な り ま す . 

この よう に プリ ンタ ・ ポー ト の 速度 が ボトル ネッ ク に な 
っ て いま す . し か し , XC3S50 の JTAG クロ ッ ク は 33MHz 
まで 上 げ る こと が で きる の で , より 高速 に TCK を 動か せる 
専用 の ハー ド ウェ ア を 作れ ば , 最大 で 100 倍 程度 の 高速 化 
が 見 込め ます . 


CO 


人 @ Pentium 4 と 比較 する と … 

690s と いう FFT の 速度 が 速い の か どう か と いう こと で 
す が , 同じ 計算 を 24GHz の Pentium 4 で 実行 し た と ころ , 
最短 で 70s 程 度 で し た . も ちろ ん 単純 な 比較 は で きま せ 
ん が , CPU の ほう が クロ ッ ク が 60 倍 速い の に 対し て 性 能 は 
10 倍 に と どまり まし た . 計算 の 効率 は FPGA の ほう が 6 倍 良 
いと いえ ます . 


@ Spartan- | と 比較 する と … 

さて , Spartan-I で 同じ 回 路 を 作っ た 場合 , どの くら い の 
性 能 に な る の で し ょ うか ? 

XC2S50 の 場合 , 16X 16 ビ ッ ト 乗算 器 符号 な し 整数 ) 1 個 
で 全体 の 約 20% の 論理 ブロ ッ ク を 消費 し ます . XC3S50 と 
XC2S50 と で は Block RAM の 構成 が 違い ます し , 乗算 の 際 
の 符号 を 考慮 し て いな い の で 正しく は 動作 し ませ ん が , 無理 
や り フ ィ ッ ト させ た と ころ , 全体 の 約 98% の 論理 ブロ ッ ク 
を 消費 し , 動作 速度 も 39MHz と な り ま し た . 

この こと か ら も , Spartan-3 の 乗算 器 が ディ ジタル 信号 処 
理 に 有効 で や る こと が わか り ま す . 


、 まとめ ) / 


実は , この ハー ド ウェ ア ・ フ ー リ 変換 機 は , 筆者 が 約 8 
年 前 に 考案 し , 標準 論理 ITC と ハード ウェ ア 乗 算 器 , GAL な 

ど , 合計 60 個 程度 の IC を 使っ て ISA バ ス の 万 能 基 板 上 に 組 
み 上 げた 複雑 な 回 路 で し 写真 2. それ が 1 チッ プ で で き 
て し まう の で すか ら 時 代 の 進歩 を 感じ ます . 

この よう に FFT を 1 チッ プ 化 する こと が で きま し た が , 
FPGA の リソー ス は まだ 44% し か 消費 し て いま せん . 窓 
数 や AD コン バー タ の イン ター フェ ー ス な ど を 除い た , 純粋 
に FEFT の 部 分 だ け な ら 30% 台 で す . 

さら に 余力 の ある 人 は , 次 の よう な 改造 を 行っ て みる と お 
も し ろ い で し ょ う . 


ET Nb 


、、 1 w て 1 
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写真 2 8 年 前 に 製作 し た ハー ドウ ェ ア ・ フ ー リ エ 変 換 機 
標準 論理 |C じ と ハー ドウ ェ ア 乗 算 器 , GAL な ど , 約 60 個 の IC を 使用 し た . 


e アン チエ イリ アシ ング ・ フ ィ ル タ の 実装 
@ 作 - シ 型 人 必 コ ン バ ー タ の 実装 
e 逆 FFT と A- 之 型 DA コン バー タ の 実装 
。oe (GT の ) の 
スペ クト ラム ・ ア ナラ イザ で 遊ぶ 際 に は , 奇声 を 発し て 周 
囲 の 迷惑 に な か ら な いよ うに , く Re だ さい . 
今回 の 設計 デー タ や 関連 ファ イル は , 本 誌 の 付属 CD-ROM 
に 収録 し て いま す . また , ソフ トウ ェ ア の アッ プ デ ー ト や 
FPGA の 回 路 の 改良 な ど は , 筆者 の ホー ムペ ー ツ http// 
wwwanahitechcom/nahitafu/ ) で 公開 し て いく 予定 で す . 
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